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Подано результати попередньої обробки первинних даних, отриманих з режимної геофізичної автоматич- 
ної станції. Як методи обробки, використано описову статистику та моделювання часових рядів. Подання резуль- 
татів багатовимірними графіками дозволило виявити феномен збігу в першому показнику описової статистики, а 
в другому - збіг коефіцієнтів моделі для подобових вимірів природного електричного поля. Використання мір 
центральної тенденції та мір варіації дає не лише основні загальні характеристики отриманої вибірки даних, але 
подання їх у вигляді багатоелементного графіка дозволяє виявити окремі специфічні точки для послідовності діб 
місяця. Їх специфічність полягає в тому, що вони вказують на доби, які за цими показниками є найбільш подібни- 
ми між собою. Використання методів описової статистики та часових рядів є доцільним для обробки первинних 
даних, оскільки їх результати визначають нові задачі для подальшого аналізу не тільки геофізичних даних, але й 
інших подібних явищ, таких як кардіологічні сигнали, показники сонячної активності та інші. 
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Вступ 

Для дослідження та вивчення фізич- 
них процесів і явищ, що відбуваються у по- 
верхневих шарах та надрах Землі, а також 
на її поверхні й в навколоземному просторі, 
використовують різноманітні геофізичні 
методи та методики. Такі дослідження ор- 
ганізовуються у різний спосіб: це можуть 


бути систематичні наукові дослідження у 
формі безперервного моніторингу чи ре- 
жимних геофізичних спостережень, або ж 
періодичні вивчення окремих регіонів чи 
площ, Кожен з таких напрямків дає важли- 
ву первинну інформацію для подальшого 
вивчення цих процесів. 
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Безперервні геофізичні спостережен- 
ня ведуться постійно в стаціонарних та 
тимчасових пунктах спостережень. На сьо- 
годнішній час, збір та реєстрація фізичних 
полів відбувається переважно автоматич- 
но, шляхом створення і поповнення різно- 
манітних баз даних. За даними регулярних 
та багаторічних спостережень можна оці- 
нювати та аналізувати параметри полів та 
їх характеристики, вивчати динаміку змін 
упродовж тривалого періоду часу, контро- 
лювати | критично-важливі | показники. 
Зазвичай, подібні дослідження забезпечу- 
ють міжнародні мережі, світові та націо- 
нальні центри, консорціуми, геологічні, 
екологічні та інші служби, об'єднання, 
агентства. Згадані організації, окрім інших 
типів даних, працюють з геомагнітними, 
сейсмічними, магнітотелуричними спосте- 
реженнями, формують відповідні архіви, 
подають візуальне представлення таких да- 
них у режимі опіїпе і, крім того, надають 
спеціальний або безпосередній доступ до 
них. 

Однак, поряд із таким масштабним 
вивченням світового рівня геофізичної ін- 
формації значний інтерес у царині геофізи- 
ки мають і локальні дослідження окремих 
регіонів чи окремих ділянок земної поверх- 
ні. Такі локальні дослідження корисні тим, 
що дозволяють глибше, точніше і детальні- 
ше вивчати конкретний регіон, досліджу- 
вати локальні зміни чи особливості, харак- 
терні саме для цієї області. Окрім цього, 
такі відомості можуть слугувати важливою 
допоміжною, додатковою інформацією для 
обробки та аналізу глобальних геофізичних 
даних. 

Постановка проблеми та формулю- 
вання мети статті 

Переважно, автоматичні стаціонарні 
пункти спостережень геофізичних даних 
представляють собою структури, утворені 
набором спеціальних сенсорів (давачів), 
реєструючої апаратури та засобів передачі 
інформації в пункти призначення. Отрима- 
ні дані кодуються відповідним чином 1 пе- 
редаються в режимі, практично, реального 
часу у відповідні організації для обробки та 
збереження. Особливістю таких даних є те, 
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що збір таких даних триває вже досить ба- 
гато часу, причому вимірювання значень 
геофізичних параметрів здійснюються по- 
стійно 1 через досить короткі проміжки ча- 
су, наприклад, у даному випадку через 
кожні 5 хв. протягом години, доби, місяця, 
року, 1 так з року в рік. Очевидно, що такі 
обсяги даних потребують спеціальних ме- 
тодів, математичного, програмного та ком- 
п'ютерного забезпечення. 

Метою дослідження є орієнтація на 
вибір простих методів первинної обробки 
даних та розробка відповідних математич- 
них моделей динаміки їх розвитку в сенсі 
створення спеціального вебсайту для від- 
повідних служб та наукових дослідників. 

Аналіз літературних джерел 

Проведений аналіз літературних дже- 
рел свідчить, що багато сайтів та організа- 
цій світового рівня і глобального масштабу 
давно працюють, аналізують і представля- 
ють результати геофізичних спостережень 
у реальному масштабі часу: Міжнародна 
мережа геомагнітних обсерваторій 
«ПКХТЕКБМАСМЕТ»  (|1,2|, Національні 
центри екологічної інформації «МСНІ», 
м. Ашвілл, Північна Кароліна, США |З3,4Ї, 
Міжнародний консорціум сейсмологічних 
інститутів «ІКІ5», США (5,6), Опіед 5аїез 
Сеоіовісаї! 5игуеу «0505» |71. Велика ува- 
га надається сейсмологічним даним як у 
нас в Україні, так і за кордоном (|5,6,98- 
12,13,14). Проте, локальні електромагнітні 
спостереження й питання обробки такого 
типу даних регіонального характеру на 
сьогодні висвітлені вкрай мало. Цікавою з 
цього приводу є робота |15|, де порушу- 
ються питання організації геоелектричного 
моніторингу в Національній сонячній 
обсерваторії зі штаб-квартирою в Боулдері 
Геологічної служби США. Прикладом ста- 
тистичного аналізу та обробки даних мере- 
жі атмосферних вимірювань є робота (16). 
Подібною за характером подання подобо- 
вих даних та їх обробки є робота (171, яка 
присвячена аналізу динаміки продажу наф- 
ти, використовуючи показники середньої 
ціни, інтервалу значень ціни, обсягу про- 
даж. Крім того, також проводиться аналіз 
та моделювання тренду подобового часо- 
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вого ряду, а візуалізація закономірностей 
представлена графічно. 

Виклад основного матеріалу 

Характерні особливості постійних 
локальних електромагнітних досліджень 
можна розглянути на прикладі РГС 
«Нижнє Селище» (Закарпатська область, 
Хустський район, географічні координати: 
48.1974729 пн.ш., 23.4560297сх д.). Станція 
розташована на території Закарпатської 
сейсмоактивної зони, тому, окрім інших 
геофізичних вимірювань, на станції 
ведуться | неперервні | електромагнітні 
спостереження. 

Для досліджень використовується 
метод природного електричного поля 
(ПЕТ), суть якого полягає у визначенні па- 
раметрів ПЕП у двох взаємно перпендику- 
лярних напрямках: північ-південь (Пн-Пд), 
захід-схід (Зх-Сх). Для цього на території 
режимної геофізичної станції розгорнута 
відповідна вимірювальна установка. Вона 


Кі - канал КІ 
ЕК» - канал К2 


Зх 


складається з двох пар електродів, що розта- 
шовані строго за напрямком Пн-Пд, Зх-Сх. 
Електроди кожної пари зафіксовані на пев- 
ній відстані, а саме: електроди напрямку 
Пин-Пд знаходяться на відстані 60 м один 
від одного, Зх-Сх - 80 м. 

Слід також зазначити, що установка є 
стаціонарною, забезпечує вимірювання 
ПЕП у двох напрямках, реєстрація інфор- 
маційних сигналів відбувається за єдиною 
методикою з використанням незмінного 
складу апаратури. Схема установки зобра- 
жена на рис. 1, принцип роботи такий: дві 
пари електродів постійно вимірюють ана- 
логові значення канальних напруг КІ, К2, 
що відповідають різниці потенціалів НІ, Е2 
на електродах напрямку Пн-Пд (ЕТ), Зх-Сх 
(Е2). У подальшому ці сигнали через 
аналого-цифрові перетворювачі перетво- 
рюються у цифрові. 


єлеЕктроди 


Рис. 1. Схема вимірювальної установки для постійних досліджень 
ПЕП на РГС «Нижнє Селище» 


Отриманий файл даних містить таку 
інформацію. В комірці А2 вказана дата 
проведення вимірів, тобто, яка доба, в ко- 
мірці В2 - назва станції, в комірках А3, ВЗ 
і С3 - специфічна технічна інформація, в 
стовпчику А, починаючи з комірки А4, по- 


даються значення годин і хвилин З 5- 
хвилинним інтервалом. У стовпчику В, 
починаючи з комірки ВА, вказується число 
місяця 19 липня Й у наступних 12 клітин- 
ках, починаючи з ВЗ, з інтервалом у 5 хви- 
лин, подаються значення електричного 
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поля з першого каналу - КІ. У комірці ВІ7 
знову вказано число місяця 19 липня й далі 
12 значень електричного поля, виміряних з 
5-хвилинним інтервалом, і так до закін- 
чення доби, причому останній відлік відпо- 
відає 23 годині 55 хвилин. У комірці СА на- 
ведене значення температури, яке слід 
читати як 23,68 "С, а починаючи з комірки 
С5, надаються значення електричного поля 
через кожні 5 хвилин з другого каналу - 
К2. У комірці СІ7 наводиться знову зна- 
чення температури й наступні 12 клітинок 


А в С в) 
1 
2 19.07.2017 М5Е| 
3 6767 
4 0-00 
5 0 к1017 
6 5 к1051 
б 10 к1082 
8 15 к1115 
9 20 к1109 
10 25 к1148 
11 30 к1119 
12 35 к1136 
13 40 к1127 
14 45 К1126 
15 50 к1113 
16 55 К1179 
Ці 1:00 
18 0 к1209 13:00 
19 5 к1204 
20 10 к1193 
21 15 К1154 
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55 К0523 


0 к0582 
5 к0588 
10 к0586 


- це значення електричного поля другого 
каналу. 

Хоча така форма є зрозумілою для 
фахівців, проте вона є досить незручною 
навіть для попередньої обробки, оскільки 
потребує додаткового перетворення в стан- 
дартні формати. Іншими словами, потрібно 
виключити неінформативні показники (го- 
дини, дати й температури) й створити для 
них відповідне подання, наприклад, у фор- 
мі окремих стовпців, певних міток тощо. 


зі (е| н І 


Тре огідіпа! Чдаїа йгот Нагам/аге 
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Рис. 2. Файл за 19.07.2017 року з оригінальними (первинними) даними природного електричного 
поля, що надійшов на сервер з режимної геофізичної станції «Нижнє Селище» 


Підготовка даних для інтерпретації 

Отримані дані мають надзвичайно ве- 
ликий обсяг, що створює значні труднощі 
для їх обробки. Лише протягом доби вибір- 
ка містить 288 варіантів, що є достатнім 
для репрезентативного статистичного до- 
слідження. Якщо розглядати дані за місяць 
(30 днів), то обсяг вибірки становитиме 
8640 значень варіантів, а за рік (365 днів) - 
103680 значень. Очевидно, аналіз таких да- 
них варто здійснювати, використовуючи 
більш короткі інтервали часу спостережень 
і загальний висновок робити на підставі 


отриманих результатів, наприклад, у межах 
доби, місяця. Тоді, зіставляючи подобові 
дані, отримуємо ситуацію протягом місяця, 
а зіставляючи щомісячні дані, отримуємо 
ситуацію протягом року 1 т.д. 

Другим моментом є те, що дані цих 
спостережень, фактично, є неперервно- 
дискретним випадковим процесом або, ін- 
шими словами, їх можна віднести до кате- 
горії нестаціонарних еквідистантних часо- 
вих рядів. У цьому випадку, обробка тако- 
го часового ряду потребує побудови мате- 
матичної моделі, яка має характеризувати 
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його загальну тенденцію та відповідні па- 
раметри його структури. "Статистичну 
експрес-обробку отриманих електромагніт- 
них даних пропонується розділити на кіль- 
ка етапів. 

Етап 1. Подання подобової динаміки 
показників описової статистики у вигляді 
загального графіка місячної тривалості. 

Етап 2. Моделювання подобової ди- 
наміки розвитку процесу з обгрунтуванням 
вибору математичної моделі. 

Етап 3. Визначення ролі динаміки ко- 
ефіцієнтів моделі, апроксимуючої добовий 
тренд, шляхом виявлення збігів моделей у 
різних точках досліджуваного процесу. 

Таким чином, у результаті цілодобо- 
вої роботи стаціонарного пункту постійних 
вимірювань електричного поля, отримують 
набір певним чином представлених даних 
електричної складової геофізичних показ- 
ників для цього місця локалізації. 

Визначення показників описової 
статистики 

Нехай отримані геофізичні дані пред- 
ставляють значення динаміки показників 
постійного електричного поля, тобто 


вибірка Х з|хі м, Хо...» Ху бо п - 288) 
характеризує його значення протягом 
доби. Найбільш поширеними в даному ви- 


падку методами попередньої обробки да- 
них є методи визначення показників описо- 
вої статистики. Ці показники описують 1 
подають основні статистичні властивості 
даних. Вони дозволяють узагальнювати 
первинні результати і використовуються 
для розв'язання таких двох задач: показати 
загальне в характері сукупності даних; по- 
казати в чому і наскільки дані різняться між 
собою. Для вирішення першої задачі нами 
використано показники центральної тен- 
денції - середнє арифметичне, медіана, мо- 
да, а для вирішення другої - показники 
варіації, в яку включено  середньо- 
квадратичне відхилення, розмах (інтервал). 
Показник дисперсії не використовувався, 
оскільки він є квадратом  середньо- 
квадратичного відхилення. Крім того, за- 
мість показника дисперсії використано 
коефіцієнт варіації. Його функція в цьому 
дослідженні - це визначення мінливості 
досліджуваного показника - значення по- 
стійного електричного поля як відношення 
середньоквадратичного відхилення до 
середнього арифметичного. 

Дані, отримані протягом однієї доби, 
мають вигляд, зображений у таблиці І. 


Таблиця І. 5-хвилинні значення природного електричного поля за добу 


Хо Години: Значення Хо Години: Значення Хо Години: Значення 
з/п хвилини на з/п хвилини на з/п хвилини на 
19.07.17 19.07.17 19.07.17 
1 0:00 110,17 21 1:40 112 я ам тт 
2 0:05 110,51 22 1:45 111,65 214 22:45 107,48 
3 0:10 110,82 23 1:50 112,42 2715 22:50 108,21 
4 0:15 111,15 24 1:55 112,57 276 22:55 108,44 
5 0:20 111,09 25 2:00 112,22 277 23:00 108,27 
6 0:25 111,48 26 2:05 112,62 278 23:05 108,61 
7 0:30 111,19 27 2:10 112,78 279 23:10 108,3 
8 0:35 111,36 28 2:15 112,96 280 23:15 108,99 
9 0:40 111,27 29 2:20 112,56 281 23:20 108,12 
10 0:45 111,26 30 2:25 112,37 282 23:25 108,84 
11 0:50 111,13 31 2:30 112,84 283 23:30 108,5 
12 0:55 111,79 32 2:35 112,87 284 23:35 108,77 
13 1:00 112,09 33 2:40 112,84 285 23:40 109,19 
14 1:05 112,04 34 2:45 113,08 286 23:45 109,13 
15 1:10 111,93 35 2:50 113,16 287 23:50 109 
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декілька добових графіків даних за поряд- 
ком їх отримання у часі. 


9 1200 1424 1648 1912 2136 00 
час (за Гриквічем) 


РГС"Н Селище" 28 07.2017р. 


ож 1200 1424 1648 1912 2136 0 
час (за 


Рис. 3. Подобова динаміка показників природного електричного поля 


Практично, форми тенденції зміни 
величини значень електричного поля є між 
собою дуже подібними, що свідчить про 
вплив невідомого на даний момент факто- 
ра, який у межах доби впливає одним і тим 
самим чином на величину електричного 
поля. Специфічною особливістю його дії є 
те, що можна досить чітко розмежувати 
його поведінку, поділивши добу на чотири, 
приблизно однакові, інтервали, тобто на 6- 
ти годинні інтервали. З графіка випливає, 


Таблиця 2. Показники описової статистики 


що в першому інтервалі значення елект- 
ричного поля монотонно зростає. У на- 
ступних двох інтервалах ці значення моно- 
тонно зменшуються, а в останньому вони 
знову зростають з тією ж самою монотон- 
ністю. У дослідженні, проведеному мето- 
дами описової статистики, отримано показ- 
ники, включаючи 1 коефіцієнт варіації для 
місяця липня 2017 року, які наведені в 
таблиці 2. 


Показники описової статистики 
Дата День Міри центральної тенденції Міри варіації 
з/п Середнє Медіана Мода Стандартне Розмах Коефіцієнт 
арифметичне відхилення | (інтервал) варіації 
01.07.2017 1 131,32 131,09 138,53 3,66 17,17 0,03 
02.07.2017 2 128,48 128,57 128,92 2,32 9,77 0,02 
03.07.2017 3 128,75 128,94 129,00 1,34 9,27 0,01 
04.07.2017 4 127,89 127,17 130,51 2,03 8,83 0,02 
05.07.2017 5 129,55 129,49 130,66 1,40 9,05 0,01 
06.07.2017 6 122,58 120,61 126,56 4,16 15,38 0,03 
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Продовження таблиці 2 


07.07.2017 7 119,63 120,30 121,29 2,25 11,61 0,02 
08.07.2017 8 120,09 120,34 120,41 2,29 17,55 0,02 
09.07.2017 9 117,87 117,40 120,03 3,27 13,37 0,03 
10.07.2017 10 114.49 115,33 116,09 3,15 15,22 0,03 
11.07.2017 11 109,67 108,68 108,45 4,61 20,57 0,04 
12.07.2017 12 112,73 113,01 113,22 1,49 9,77 0,01 
13.07.2017 13 112,30 112,39 112,17 1,83 14,58 0,02 
14.07.2017 14 110,09 109,82 114,06 3,09 12,31 0,03 
15.07.2017 15 110,46 110,78 111,42 1,54 8,27 0,01 
16.07.2017 16 108,40 108,93 111,60 2,59 12,09 0,02 
17.07.2017 17 107,73 108,11 111,15 2,94 15,55 0,03 
18.07.2017 18 107,96 108,66 110,23 2,18 22,89 0,03 
19.07.2017 19 108,28 108,03 112,96 3,57 16,18 0,03 
20.07.2017 20 103,57 102,37 111,53 5,96 25,62 0,05 
21.07.2017 21 100,76 99,62 98,66 5,11 18,96 0,05 
22.07.2017 22 94,04 92,56 99,76 6,03 27,07 0,06 
23.07.2017 23 88,01 86,45 96,95 6,91 20,58 0,08 
24.07.2017 24 79,41 77,05 87,89 7,48 45,59 0,09 
25.07.2017 25 75,01 75,63 69,92 3,40 11,28 0,05 
26.07.2017 26 75,87 76,11 76,35 2,46 10,81 0,03 
27.07.2017 27 77,10 79,02 81,06 3,78 12,16 0,05 
28.07.2017 28 70,54 71,87 75,16 431 13,35 0,06 
29.07.2017 29 68,53 68,40 68,24 2,59 9,54 0,04 
30.07.2017 30 65,88 66,38 71,88 5,92 18,26 0,09 
31.07.2017 31 65,26 65,56 65,79 4,78 16,93 0,07 


Отримані результати характеризу- 
ють центральну тенденцію та варіацію да- 
них протягом одного місяця. Значення ін- 
ших показників описової статистики, та- 
ких як стандартна помилка, асиметрія, 
ексцес, мінімум, максимум, сума та кіль- 
кість у цьому дослідженні не використо- 
вувались. Справа в тому, що стандартна 
помилка тут визначається відношенням 


с 
п 2 
відхилення та кількість даних відповідно. 
Оскільки для всіх подобових даних кіль- 
кість однакова, то в зазначеному випадку 
цілком достатнім є показник варіації 
середньоквадратичного відхилення. 
Показники асиметрії та ексцесу по- 
в'язані з функцією закону розподілу. 
Оскільки в цьому дослідженні не стояла 
задача визначити закон розподілу цих да- 
них, тому й не було підстав уточнювати ці 


де с, п - середньоквадратичне 


показники. Щодо показників максимуму й 
мінімуму, то вони на зазначеному етапі до- 
сліджень вказують на присутність екстре- 
мальних значень. Показник суми - суму 
значень варіантів за однакової кількості 
варіантів для кожної доби можна потракту- 
вати як добову потужність сигналу. Показ- 
ник кількості даних відомий і однаковий 
для кожної доби. 

На рис. 4 зображено подобові графі- 
ки. На графіках дуже чітко виражена сут- 
тєво нелінійна зміна показників централь- 
ної тенденції. Для неї властивою є точка 
перегину, що означає присутність моди, 
що й підтверджують дані описової статис- 
тики. Важливим візуальним моментом є 
те, що характер зміни варіантів дуже 
подібний. 
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Динаміка показників описової статистики за місяць липень 2017 р 


12 3 4 5 6 7 8 9 101112 13 14 15 16 17 18 1920 21 22 23 24 25 26 27 28 29 30 31 
Подобовізначення 


Рис. 4. Динаміка показників описової статистики протягом одного місяця 


За результатами подобових показни- 
ків описової статистики побудовано бага- 
товимірне графічне їх подання. Врахову- 
ючи суттєву відмінність значень показни- 
ків середнього, моди, медіани від показни- 
ків  середньоквадратичного відхилення, 
розмаху та коефіцієнта варіації, усі зна- 
чення були приведені до одиничного ін- 
тервалу у своїх вибірках за формулою: 
оо ті | де Хє Го, 1 - приведене 

Хтах 7 Хтіп 
до одиничного інтервалу значення х;,. 


Х 


Результати такого нормування подані у 
вигляді багатовимірного графіка на 
рис. 4. 

Подання результатів описової ста- 
тистики в такий спосіб, як це зроблено на 
рис. 4, не є традиційним, проте його вико- 
ристання виявило певні збіги показників 
описової статистики для деяких діб. На- 
приклад, для показників центральної тен- 
денції особливими можна вважати доби: 
1, 6, 14, 19, 29, хоча більш точно їх можна 
визначити за табл. 2. Для показників варі- 
ації, такими є 3, 5, 12, 15, 29 доба. Пове- 
дінку рівнів (абсолютних значень вимірів 
електричного поля) протягом місяця мож- 
на вважати «квазісинхронною», оскільки 
форма тренду для кожної доби фактично 


У()4у у Зо)» 1 


впорядкованих у часі випадкових величин 
і характеризувати (визначити, описати) 


його як лінійно  нестаціонарний, з 
добовою періодичністю, тривалістю 
100 


| ер 


повторюється, проте виняток може вклю- 
чати лише одну-три доби. 

Існування таких збігів може бути ви- 
падковою подією, але може бути й резуль- 
татом відсутності будь-яких впливів, тобто 
можна вважати стан середовища нормаль- 
ним, а усі інші - як результат впливу різ- 
них факторів. Варіантів їх трактування 
може бути дуже багато, а тому досліджен- 
ня в цьому напрямку мають певний 
науковий сенс для отримання наукового 
результату. 

Присутність таких особливостей у 
даних вимагає нових досліджень даної 
місцевості та врахування факторів, які за- 
відомо можуть впливати на значення при- 
родного електричного поля, наприклад: 
атмосферні, кліматичні, гравітаційні, сейс- 
мічні, магнітні тощо. 

Аналіз трендів подобових часових 
рядів 

Отримані дані можуть бути пред- 
ставлені у вигляді графіка еквідистантно- 
го часового ряду, зображеного на рис. 5. 
Графічне зображення дає загальний ви- 
гляд його структури, вказує на характер 1 
вид стаціонарності та коливальну складо- 
ву. Цей ряд можна формально подати як 
множину 


288, 1; - бтіп, з єТ) 


ДІ - 288: 5тіп - 24йоийг, де Т - множина 


5-хвилинних інтервалів часу. На рис. 5 
наведені графіки для діб липня 2017 року. 
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Значення природного електричного поля на липень 2017 р. 
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Рис. 5. Значення постійного електричного поля протягом липня 2017 року 


Оскільки форма рівнів часового ряду 
повторюється, відмінності її поведінки 
можна оцінити, використовуючи значення 
коефіцієнтів моделі, якою апроксимують 


01.07.2017 


дтричав 1, тм 


12.07.2017 


дтрічів ЕТ, тм 


дтрімав Е1, т 


дтрійцав Е1, т 


добові рівні. У цьому випадку, візуально, 
тренд добових значень має точку переги- 
ну, яка переважно припадає на середину 
доби, як відображено на рис. б. 


1912 


23.07.2017 


Рис. 6. Тренди подобових даних, апроксимовані поліномами третього ступеня 


На перший погляд, такою апроксиму- 
ючою функцією у(г) в якості моделі 
тренду доцільно випробувати поліном 
третього ступеня: 


у(г)е аз б аг 4 ар ад. 
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Саме властивість поліному третього 
ступеня мати можливість відтворювати 
залежність, яка має точку перегину, робить 
його найпростішою функцією для апрокси- 
мації цих даних. Основним недоліком полі- 
номів є те, що практично неможливо надати 
їх коефіцієнтам смислового навантаження 
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чи фізичного тлумачення. Тим не менше, дані значення коефіцієнтів апроксимуючих 
зроблено спробу такої апроксимації. На поліномів. 
рис.6 наведені графіки апроксимованих 
трендів подобових даних. У таблиці 3 по- 
Таблиця 3. Результати апроксимації даних для кожної доби 


Доби з/п Коефіцієнти поліному к2 Середнє У тодеї 
аз а2 а! ад 
1 112,83 -160,06 50,399 131,27 0,7843 131,32 131,3193 
2 41,283 -52,377 9,3432 130,94 0,7488 128,45 128,4771 
3 15,403 -13,579 0,0221 129,42 0,3302 128,75 128,7533 
4 40,329 -44,379 4,3725 130,41 0,7747 127,89 127,8959 
5 30,727 -38,171 9,4932 129,84 0,3885 129,55 129,5422 
6 62,384 -76,276 7,148 128,52 0,9198 122,58 122,5756 
7 63,741 -79,745 18,979 120,79 0,7487 119,63 119,6289 
8 17,52 -7,2516 -7,1141 121,69 0,4313 120,09 120,0914 
9 33,85 -39,632 -0,8752 123,05 0,8976 117,87 117,8758 
10 64,724 -91,144 25,436 115,97 0,6381 114,49 114,4897 
11 70,745 -75,838 0,4634 117,02 0,8674 109,67 109,6669 
12 46,178 -69,773 26,774 111,06 0,4236 112,73 112,7288 
13 -11,221 31,493 -22,318 115,77 0,3997 112,30 112,3073 
14 64,952 -79,563 13,084 113,83 0,9097 110,09 110,092 
15 34,995 -40,494 8,162 111,13 0,5412 110,46 110,4578 
16 45,503 -49,228 3,8593 111,5 0,6568 108,40 108,3963 
17 78,904 -105,7 27,344 109,57 0,8803 107,73 107,734 
18 90,163 -123,94 37,578 107,92 0,7505 107,96 107,9636 
19 96,442 -130,49 34,345 110,49 0,9241 108,28 108,2761 
20 94,054 -102,42 6,0453 111,17 0,7465 103,57 103,5705 
21 114,98 -158,75 -158,75 104,72 0,9284 100,76 100,7684 
22 124,53 -163,79 35,526 99,73 0,8388 94,04 94,03498 
23 178,32 -246,57 66,487 92,37 0,9637 88,01 88,00613 
24 185,08 -267,01 80,343 81,971 0,8334 79,41 79,41108 
25 32,683 -55,933 15,123 77,175 0,924 75,01 75,00806 
26 46,827 -74,739 39,361 69,163 0,8854 75,87 75,36668 
27 78,867 -122,61 40,22 78,137 0,9224 77,10 77,69869 
28 100,35 -144,68 41,727 72,815 0,9246 70,54 70,5427 
29 96,522 -141 48,317 67,248 0,9305 68,53 68,52966 
30 180,7 -259,95 81,759 66,446 0,925 65,98 65,37847 
31 178,78 -261,3 90,269 62,712 0,9228 65,26 65,26173 
Як і у випадку з показниками описо- тут також існують специфічні точки, в яких 
вої статистики, був використаний підхід -- абсолютні значення коефіцієнтів а:;, де 


подання значень коефіцієнтів поліному на 
одному графіку, зображеному на рис. 7. 
Особливість результатів полягає в тому, що 


7 - 0, 1, 2, 3, мають дуже близькі значення. 
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Подобова динаміка значень коефіцієнтів моделей 


Коефіцієнти моделі 


т т т п Т т т ї 
б 171:18|19:20:21|221:231|24 


Номер доби 


Рис. 7. Зміна коефіцієнтів моделі щодоби 


Оскільки другий коефіцієнт а» є 
від'ємним, у результаті цього спостеріга- 
ється симетричність графіків для коефіці- 
єнтів а» і аз. Якщо подати цей графік не з 
реальними значеннями, а з абсолютними, 
тобто представити значення коефіцієнтів 
, трафік 
матиме вигляд, зображений на рис. 8. 

Вже візуальний аналіз графіків на 
рис. 8 вказує на особливі точки, в яких 
значення коефіцієнтів моделі збігаються. 
Такі збіги є нерегулярними й неабсолют- 
ними, тобто повний збіг не завжди має 
місце. Для таких діб як 3, 8, 13 125 збіги є 
досить близькими. Природу цих збігів 


а; як абсолютні значення а і 


важко пояснити, використовуючи лише ці 
дані. Очевидно, що такі збіги, в першу 
чергу, можна вважати випадковими (як 
збіг випадкових обставину). 

Коефіцієнт ад є дуже близьким до 


середнього арифметичного для даної доби. 
Иого поведінка, як видно з графіка на 
рис. 7, є досить монотонною і відносно 
гладкою. Порівняно з іншими місяцями, 
прослідковується певна  коливальність 
цього показника з періодом у декілька мі- 
сяців. Цей феномен також потребує відпо- 
відних досліджень. 


Подобова динаміка абсолютних значень коефіцієнтів 
моделей 
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Рис. 8. Динаміка абсолютних значень коефіцієнтів моделей 


Подання поведінки коефіцієнтів а; 
моделі на одному графіку вказує на те, що 
аз 1 аз змінюються пропорційно, 1 на дея- 
ких етапах їх повторює показник а . Отже, 


для деяких діб існує подібність у структу- 
рах апроксимуючих моделей. Оскільки 


якість апроксимації за коефіцієнтом детер- 


мінації к2 0.89 , можна вважати, що від- 
повідність моделі добовим значенням є до- 
статньою для використання цього підходу 
-- встановлення специфічних точок динамі- 
ки природного електричного поля для да- 
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ної локальної режимної геофізичної авто- 

матичної станції «Нижнє Селище». 

Висновки 

Проведення попередньої обробки 
первинних геофізичних даних, автоматич- 
но переданих з режимної геофізичної стан- 
ції щодобовими специфічними файлами, 
підтвердило доцільність використання 
таких підходів. 

1. Використання мір центральної тенденції 
та мір варіації дають не лише основні 
загальні характеристики отриманої ви- 
бірки даних, але їх подання у вигляді 
багатоелементного графіка дозволяє ви- 
явити окремі специфічні точки для по- 
слідовності діб місяця. Їх специфічність 
полягає в тому, що вони вказують на 
доби, які за цими показниками є най- 
більш подібними між собою. Цей фено- 
мен можна трактувати по-різному: від- 
сутність зовнішніх впливів, у ці доби 
відбувся один і той самий вплив і тому 
подібне. Цілком зрозуміло, що на усьо- 
му проміжку отриманих за весь час ро- 
боти станції даних, цей феномен мусить 
мати якесь пояснення в майбутніх 
дослідженнях. 

2. Побудова подобового часового ряду з 
отриманих даних показала, що дані в ін- 
тервалі доби мають дуже часто повто- 
рюваний за формою тренд. Для аналізу 
цього явища подобові рівні часового ря- 
ду були апроксимовані поліномом тре- 
тього ступеня. Багатовимірний часовий 
ряд, утворений з коефіцієнтів цього 
поліному-моделі також показав збіги їх 
значень для деяких діб. Ці збіги, хоча і 
не відповідають таким для описової ста- 
тистики, проте нехтувати ними, не зна- 
ючи причин, вважаємо, не варто, оскіль- 
ки вони несуть якусь інформацію, мож- 
ливо, дуже важливу. 

Таким чином, можна стверджувати, 
що використання методів описової статис- 
тики та часових рядів є доцільним для об- 
робки первинних даних, оскільки їх ре- 
зультати визначають нові задачі для по- 
дальшого аналізу не тільки геофізичних да- 
них, але й інших подібних явищ, таких як 
кардіологічні сигнали, показники сонячної 
активності та інші. 
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Р.М. Камінський, Н.Б. Шаховська, 
Л.Г. Савків, Я.Ю. Варецький, 
С.В. Саварин 

Ргейтіпагу ргосеззіпя ої ргітагу єе0- 
рпубісаї дака, ацопайсаПу (гап5тійед бот 
Фе гесіте яеорпузісаї згайоп (о Ше даїйу 
зресійс Піе5 сопіттеай Фе Геаз1ір Шу ої и5іпе 
зисп арргоаспез. Тре ригро5е ої Ше гезеагсП 
15 їо Їоси8 оп Фе сроїсе ої 5ітріе птеШод5 ої 
ргітагу Чаіа ргосез5іп8, апа Ше деуеіортепі 
ої арргоргіаїе плафетайсаї птодеїз ої Бе ду- 
патіс5. Тре гезиіїз ої ргейтіпагу ргосез5іп8, 
ої Фе ргітагу аа обіаїпед гот Ше гебіте 
георпузіса! ашіотайс 5іайоп аге 5ибплійеа. 
Гезсгірйуе згайяйся апа 5ітиіайоп ої те 
5егіе5 ууеге ц5ед а5 ргосе55іпея, плеїподз8. Тре 
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ргезепіайоп ої (ре гез5иії5 бу плиїдїтепзі- 
опа! єгарр5 аПомуеа іо геусеаї Фе ррепотепоп 
ої соіпсідепсе іп (ре Пг5: іпдісаїого ої Че- 
5сгірйуе зіайяйся, апа їп Ше 5есопа, Ше соїп- 
сіЧепсе ої Ше піоайвідгароє;5 соеНйсіепія Їог 
Фе Чайу плазигетепіз ої (Фе пашга! еіесітіс 
Пе!4. О5іпє Фе теазигез ої Ше сепіта! (гепа 
апа Ше Чеєтее ої уапайоп 51уе5 пої опіІу Ше 
Базіс єепегаї сПагасіегізйїс5 ої Ше гесеїмед 
затріе ої Чака, биї Феїг гергезепіайоп іп Фе 
Гога. ої а паціб-еісетепі єгарп аПом/8 удси (о 
ідепіїу зресійс роїіпі5 їог Ше 5едиепсе ої 
Чауз ої Фе пттопіЮ. ТРеїг зресійсіку Пез іп Фе 
Гасі ШФас Феу роїпі (о Ше Чау5 Ша аге Ше 
тові 5іптіаг атопе ШПе5е іпдаїсакогя. ТРі8 рре- 
потепоп сап Бе іпіегргеїсда дійегепйу: Ше 
аб5епсе ої ехіегпа іпПиепсе5, іп ШФоз5е Чауб 
Феге у/а5 опе апа Ше 5ате іпПиепсе апа Ше 
ПКе. П 15 сіеаг Шаг різ ррепотепоп ти5і 
рауе 5оте ехріапайоп їп Гиіиге 5їшаїе5 оуег 
Фе епіге бар обіаїпед дигіпє Ше епіїге м/огк 
йте ої Ше Чака 5кайоп. Соп5ігисйоп ої Фе 
Чайу Це 5егіе5 гоп їбе Чаїа обіаїпед 
5похуса Шаг Фе Чака іп Ше іпіегуа! ої йте 
раме уегу ойеп гереййуєе Їогтп ої Ше (гепа. 
Бог Фе апаіузіз ої (різ рпепотепоп, їФе Іе- 
уеі5 ої Ше те 5егіе5 меге арргохітагеа Бу 
Фе роіупотіа! ої (пе ШФігд деєгее. А тиій- 
фітепзіопа! те 5егіе5, Гогтед фопі Ше 
соеййсіепі5 ої із роіупопиаї піодсі, аізо 
5помуеа Ше соїпсідепсе ої ШФеїг уаЇце5 Їог 5е- 
уегаї Чдауз. Тре5е соіпсідепсе5 ШоцеЮ дО пої 
сотезропа (о 5исП Фезсгірйуе зсацяйся, Биї 
пеєіесі Шет уліпоці Кпоууїпе, Ше геа5оп5 ууе 
до пої фріпК, ресац5е Шеу сапу 50те іпіог- 
тайоп, плаубе уегу штрогіапі. ТПиз, 1ї сап Бе 
агоцед Шаї Ше и5е ої де5сгірйуе 5(ай5йся апа 
йте 5егіе5 теїпод5 15 арргоргіаїе Їог Ше рго- 
сеззіпє ої ргіплагу Чака, 5іпсе ШПпеїг ге5иії5 
депе пеуу іа5К5 Їог Гигірег апаіузі5 ої пої 
опіу беорпузісаї Чака Бик аї5о оїбег 5ітийаг 
рбпепотепа 5исп а58 сагіас 5ієпаї8, 50Їаг 
аспуту апа оег5 
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